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Abstract 

We report the experimental findings of formation and motion of heap in granular materials in an 
inclined and vertically vibrated container. We show experimentally how the transport velocity of 
heap up container is related to the driving acceleration as well as the driving frequency of exciter. 
An analogous experiment was performed with a heap-shaped Plexiglas block. We propose that 
cohesion force resulted from pressure gradient in ambient gas plays a crucial role in enhancing 
and maintaining a heap, and ratchet effect causes the movement of the heap. An equation which 
governs the transport velocity of the heap is presented. 

PACS numbers: 45.70.Mg, 45.70.-n, 47.20.-k 



1 



It is well known that under vibration many processes, e. g. segregation, convection, heap- 
ing, density wave, and anomalous sound propagation, which govern the physics of granular 
materials, are quite unusual 111, so that the properties of such materials are not well under- 
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stood. For example, heaping is one of long-standing problem since Faraday [Si y, 1^ |5|] , and 
several physical mechanisms have been identified as possible causes of it: friction between 
the walls and particles , analog of acoustic streaming if the shaking is nonuniform , gas 



pressure effect j8[, and auto-amplification [9||. Recently we have observed experimentally 
the formation of a heap and the motion of it from lower to higher end of an inclined and 
vertically vibrated container. The experiment was conducted in a Plexiglas rectangular con- 
tainer [370 mm (length) x25 mm (width) x80 mm (height)]. We investigated the behavior 
of two types of quartz sands: spheres of diameter 0.15-0.20 mm and grains of irregular shape 
or coarse surface with diameter 0.3-0.5 mm. The container was inclined with an inclination 
a from 0.04 radian to 0.25 radian by putting a pad underneath it. The vibration exciter 
(Briiel & Kjaer 4805) was driven by a sinusoidal signal, and controlled by a vibration exciter 
control (Briiel & Kj^r 1050). Driving frequency / and dimensionless acceleration amplitude 
r = Aix'^pA/ g (where A is driving amplitude and g the gravitational acceleration) were used 
as two control parameters. 

The experiment shows that even if in horizontal container the horizontal acceleration 
can cause movement of a heap. The horizontal acceleration of our exciter is about 2.5% 
of the vertical acceleration. To get rid of the influence of the horizontal component of the 
driving acceleration on the movement of the heap along the longer direction of the container, 
we adjusted the longer direction of the container orthogonal to horizontal acceleration of 
the exciter. The ranges of F and / we used were from 1.4 to 2.8 and from 11 Hz to 20 
Hz respectively, which are good ranges for heap formation. At first, about 80 ml of sands 
was uniformly put into the lower part of the vessel. Then as F increased to and beyond 




FIG. 1: The center-high heap (a) and the wall-high heap (b) of coarse sands, a = 0.087 radian, 
/ = 15 Hz, r = 2. 
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FIG. 2: The transport velocity V of heap formed by coarse sands, (a) ^ vs F for different 
frequencies. The soHd Unes are fits hy V = A{T — Tc)^a'^ + B with a = 0.045 radian, (b) 
{V — B)/{Aa'^) vs F — Fc for the data in (a). The solid line is (F — Fc)^. (c) V vs a for different 
frequencies. The sohd hnes are fits hy V = A{T - Fc)Pa'' + B with F = 2. (d) {V - B) /[A{T - Fc)^] 
vs a for the data in (c). The solid line is a'^. (e) {V — B)/A vs (F — Fc)^a'^ for the data both in 
(a) and (c). The solid line is (F — Tf^^a'^. 

some critical acceleration Tc{> 1), center-high heaps formed, meanwhile they moved up 
the container. As center-high heaps reached higher end wall of the container, they moved 
forward continuously until they became wall-high heaps. Fig. 1 shows the photos of heaps 
formed by coarse sands. Fig. 1(a) is center-high heap, and (b) is wall-high heap. For center- 
high heap, the back (or right shown as in the figure) surface is longer than the frontal (or 
left) one, but the dynamical angles of repose (slightly smaller than the maximum angle of 
repose of the static heap) of both frontal and back surfaces relative to horizontal are the 
same. The difference between the heaps formed by two types of sands is: the dynamical 
angle of repose of the heap formed by coarse sands is larger than that formed by spherical 
sands. The heaps moved up the container with nearly uniform velocity. We measured the 
velocities V of the heaps. Fig. 2 gives the results for coarse sands. Fig. 2(a) shows ^ vs F 
for different /. One can see that V increases with F for all of frequencies, but decreases as 
/ increases for all values of F. The velocity of heap formed by coarse sands is greater than 
that of heap formed by spherical sands for all sets of values of F and /. Fig. 2(c) shows the 
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velocity of heap as a function of inclination of the container. It is shown that V increases 
as a, and behaves in the same way as in Fig. 2(a) as / changes. The heap formed by coarse 
sands moved faster than the heap formed by spherical sands did for the same inclination 
and the same driving parameters. 

Obviously, the transport of the heap is a cooperative behavior of the granular materials. 
To examine this idea we put a Plexiglas block the same in shape and dimension as the 
sand heap on the same vibrated inclined container, a similar transport of the block up the 
container was observed. So we consider that the transport of the heap is similar to that of 
a sohd block. But why does it move, or what is the mechanism of the transport? We used 
a high speed camera (Redlake MASD MotionScopc PCI 2000 SC) to record the movement 
of the block as it moved up the container with record rate of 250 fps (frames per second), 
then played back slowly (25 fps). In this way, we can see the detail of the movement of 
the block. We marked the center of mass of the block with a black point C (shown as in 
Fig. 4). Fig. 3(a) is the orbit of the center of mass of the block measured experimentally 
in laboratory reference frame. Fig. 3(b) is a schematic diagram of Fig. 3(a), with which 
we can discuss the movement of the block conveniently. To examine the effect of ambient 
gas, wc performed two things. First, to get rid of the pressure difference between the gap 
and the above of the block, we drilled a number of holes vertically and parallel each other 
through the block, and observed no transport of the block. And second, we evacuated the 
air from the container and also observed no transport of the block (either with or without 
holes drilled through the block) . These imply that air pressure difference plays a critical role 
for the movement of the block up the container. Therefore we propose a mechanism for the 




FIG. 3: (a) The orbit of the center of mass of the block measured experimentally, a = 0.06 radian, 
/ = 15 Hz, r = 2.3. (b) The schematic diagram of (a). 
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movement of the block up the container as follows. In each cycle, when T cos27r ft < —g, 
block separates with container, and a gap forms between block bottom and the container 
floor. The pressure pi in gap is less than the atmospheric pressure Pq above the block on 
an average (ref. Q] shows that the mean pressure in a gap between the bottom of granular 
bed and the container floor is below atmospheric pressure. And as Faraday Said: in this 
gap, "it forms a partial vacuum". We consider this also suits to the case of block). This 
pressure difference causes a force exerting on the block. Fig. 4 is a schematic diagram for 
the forces acting on the block (or heap) during free flight. The forces acting on the frontal 
(I) and back (II) parts (with dashed line as intersection) are represented by vectors Fi and 
F2 at the centers of mass of two parts, Ci and C2, respectively. They are perpendicular to 
frontal and back surface (i. e. the frontal and back sides in the figure), respectively. The 
magnitudes of them are proportional to lengths of the frontal and back sides, respectively. 
They are decomposed into two components: those parallel (represented by Fic and F2c) and 
those perpendicular (represented by Fid and F2d) to the bottom of the block, respectively. 
The parallel components Fic and F2C equal in magnitude but opposite in direction. So 
the total force F is only the sum of perpendicular components Fid and F2d, acts at the 
center of mass of the block, C, is perpendicular to the bottom of block, and makes an angle 
a (i. e. the inclination of the container) with gravitational force of the block. This force, 
together with gravitation force Mg (M is mass of the block), force the block to move along 
a ballistic trajectory until collides with the container [A -—>■ B in Fig. 3(b)]. Upon colliding 
with the container, the block slides down the container [B ^ C in Fig. 3(b)] due to gravity, 
meanwhile moves together with the container in vertical direction until next separation with 
container [C ^ D in Fig. 3(b)]. Then a new cycle begins. But due to the friction between 
the block and the container floor, the distance up the container during free flight is much 




atmospheric pressure 



FIG. 4: A schematic diagram for the forces acting on block (heap) during free flight. 
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larger than the distance down the container during bed-floor collision, i. e. the block moves 
one step up the container in each cycle. The analysis above shows that the transport of 
block up the container is a ratchet effect caused by the pressure difference (between the gap 
and above of the block) and the friction force (between block and the container floor). 

For the granular bed, if for any reason, some small (or flatter) initial heap has formed. In 
the period of free flight, the pressure difference between the gap and above of the heap [l^ 
leads to a pressure gradient (or force) in the interior of the heap. Also cf. Fig. 4, and we 
still use Fi and F2 as representatives of the total force acting on frontal and back parts of 
the heap, respectively. These two forces are perpendicular to the frontal and back surface, 
respectively. Here we call the parallel components Fic and F2C as a couple of cohesion force, 
i. e., they make granular materials cohesive and enhance the heap. Under the action of the 
cohesion force, heap is compressed in direction parallel to heap bottom while elongated in 
direction perpendicular to heap bottom. And this makes the heap bottom convex, as we 
have observed experimentally in both inclined and horizontal containers P|. Fig. 5(a) is a 
schematic diagram of the velocity (or mass flow) field in the interior of the heap formed 
by coarse sands during free flight. The total force F acting on the center of mass of heap, 
together with gravitational force Mg, force the heap as a whole to move along a ballistic 
trajectory and up the container. Upon colliding with the container, the center part of the 
heap bottom touches floor first, then the other parts, from center to outer, touch floor 
consecutively. This results in a further enhancement of the heap. In this way the slope of 
the heap is getting larger and larger. Once the slope angle reaches and exceeds the dynamic 
angle of repose of the heap, and upon colliding with container, the avalanche occurs on 
the surface of the heap. Then the process repeats periodically: during the free flight, the 
compression due to cohesion force makes the slope of heap exceeding the dynamic angle of 
repose, and during the bed-floor collision, avalanche occurs. So when the heap reaches a 
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FIG. 5: The schematic diagram of the velocity (or mass flow) field in a heap formed by coarse 
sands during the period of free flight (a) and heap- floor collision(b), respectively. 
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steady state, in the laboratory reference frame, one can see a steady convection flow: in the 
interior of the heap, grains move upward, while at the surface grains move rapidly downward 
[schematically shown as in Fig. 5(b)]. This analysis is also suitable to wall-high heap, in 
which the cohesion force points higher end wall. Similar to the block, the transport of the 
heap up the container is also a ratchet effect, which is caused by the pressure gradient in 
the heap and the friction force between the heap and the container floor. If we pump the air 
out of the container, as pressure is reduced, heap reaches a flat state gradually. This shows 
that ambient gas plays a important role in enhancing and maintaining a heap. 

Let us now describe the ratchet effect little more concretely but qualitatively [also cf. 
Fig. 3(b) and Fig. 4]. We denote (3 as the ratio of free-flight time to the excitation period 
T in each cycle. In each cycle, the up-container distance of free flight by heap is Si ~ 
ai/3^T^/ cos a, where a is the inclination of the container, and ai is the mean acceleration 
of the center of mass of the heap due to the horizontal component Fx of the force F. Then 
ai ~ F sma {F = |F|), and si ~ F/5^T^tana. In period of bed-floor collision, heap slides 
down the container. The distance of sliding is S2 ~ 02(1 — jS^T'^, where a2 is the mean 
acceleration down the container, which is determined by the friction coefficient fi between 
the granular bed and the container floor, the inclination a of the container, and gravitational 
acceleration g through 02 = (7 sin a — fig cos a. The total (or net) displacement of heap up 
the container in each cycle is s = Si — S2, timing driving frequency / gives the velocity V of 
heap up the container. The ratio P increases with driving amplitude A = Tg/Aix'^p i. e. 
V increases as F, but decreases as / increases. The larger the /x, the smaller the S2, and then 
the lager the V . Because the friction between coarse sands and the container floor is larger 
than that between spherical sands and the container floor, the transport velocity of heap 
formed by coarse sands is larger than that of heap formed by spherical sands. In view of 
S2 « si and a « 1 in our experiment, then V increases as inclination a of the container. 
The pressure gradient or force in the interior of heap depends on driving acceleration and 
frequency, inclination of the container, and properties of grains in a specific and complex 
way, and changes with time and space. The investigation on this problem is under way. 

Any way, from the description above we have already had a qualitative picture on how 
the transport velocity of the heap up the container depends on driving acceleration and fre- 
quency, inclination of the container, properties of the grains, and friction between heap and 
container floor. Based on this and on the analysis of experimental results, the experimental 
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data in Fig. 2(a) and (c) are well fitted by the equation 

v^A{r- r^Ya'^ + B, (1) 

with p = 0.6, q = 0.8, and a = 0.045 radian for Fig 2(a) and F = 2 for Fig 2(c), respectively. 
The solid lines in the figure are nonlinear least-square fits by Eq. (1). The parameters A, B 
and Fc are depend on the driving frequency /, the dynamic angle of repose of heap which 
depends on properties (e. g. size and shape etc.) of grains, and the friction coefficient // 
between granular bed and container floor. If the data of Fig. 2(a) are plotted as {V — 
B) / {Aa'^) vs F — Fc, they collapse onto a single curve. Fig. 2(b); and if the data of Fig. 2(c) 
are plotted as {V — B)/[A{T — Fc)^] vs a, they also collapse onto a single curve. Fig. 2(d), 
within the experimental resolution. Moreover, if the data of both Fig. 2(a) and Fig. 2(c) are 
plotted as (y — B)/A vs (F — Fc)^q;^, they all collapse onto a single curve. Fig. 2(e). 

This model is also suitable for the heap (either the center-high or the wall-high heap) 
formed in a horizontal container, where the cohesion force (as above) plays the same role 
(enhancing and maintaining the heap) as in an inclined container. The total force F due to 
the pressure gradient is perpendicular to heap bottom, i. e. is parallel to gravitational force 
Mg, and no force forces heap to move in horizontal direction. 

Our conclusion is: pressure gradient in ambient gas, which makes cohesionless granular 
materials cohesive, plays a crucial role in enhancing and maintaining a heap in vibrating 
granular materials; pressure gradient in ambient gas and friction force between the granular 
bed and the container floor, which lead to a ratchet effect, are unique cause for the transport 
of a heap up an inclined container. The transport velocity of the heap is well described by 
the equation (1). Our mechanism for enhancing and maintaining heap also shows that one 
reason for lacking of heaping in molecular dynamics (MD) simulation may be that those 
models did not take into account the cohesion force due to ambient gas effect. 
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